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LaRC IR Retrieval Algorithm 


PART A: REGRESSION RETRIEVAL (Zhou et al., GRL 2005) 

Using an all-seasonal-global training database to diagnose 0-2 cloud layers from training relative 
humidity profile: 

A single cloud layer is inserted into the input training profile. Approximate lower level cloud using 
opaque cloud representation. 

Use parameterization of balloon and aircraft cloud microphysical data base to specify cloud effective 
particle diameter and cloud optical depth: 

Different cloud microphysical properties are simulated for same training profile using random 
number generator to specify visible cloud optical depth within a reasonable range. Different habitats 
can be specified (Hexagonal columns assumed here). 

Use LBLRTM/DISORT “lookup table” to specify cloud radiative properties: 

Spectral transmittance and reflectance for ice and liquid clouds interpolated from multi-dimensional 
look-up table based on DISORT multiple scattering calculations. 

Compute EOFs and Regressions from clear, cloudy, and mixed radiance data base: 

Regress cloud, surface properties & atmospheric profile parameters against radiance EOFs. 

PART B: 1-D VAR. PHYSICAL RETRIEVAL (Zhou et al., JAS 2007) 

A one-dimensional (1-d) variational solution with the regularization algorithm (i.e., the minimum 
information method) is chosen for physical retrieval methodology which uses the regression 
solution as the initial guess. 

Cloud optical/microphysical parameters, namely effective particle diameter and visible optical 
thickness, are further refined with the radiances observed within the 10.4 pm to 12.5 pm window 
region. 
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LaRC Algorithm Flowchart 
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Hybrid Retrieval Algorithm Flowchart 
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cloud heights, and (4) water cloud 



Historical training data 
with cloud parameters for 
radiance simulation including 
a realistic cloud radiative 
transfer calculation 


Regression retrieval using “all 
cloud” coefficients category to 
predetermine cloud height 


Regression retrieval using “cloud 
height grouped” coefficients for ice 
to iterate and finalize the retrieval 


Regression retrieval using “cloud 
height grouped” coefficients for 
water to iterate and finalize the 
retrieval 



Regression retrieval using “clear 
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Cloud optical thickness and particle size 
are refined with cloud effective micro- 
window radiances: 11.1 pm to retrieve 
cloud optical depth, and 11.1 pm & 12.5 
pm windows to retrieve cloud particle 
size. 


Compute Jacobian Matrix: analytical 
Jacobian for thermodynamic parameters; 
and numerical Jacobian for cloud 
parameters (i.e., t and P cld ). 

Matrix Inversion Solution: the 
regularization algorithm (or minimum 
information method) in a nonlinear least- 
squares spectral radiance fitting mode. 




Statistical EOF Physical Regressions 




1 = 1+1 


General Matrix Inversion Solution: 

5R = A5Q, where A- = dR/dQj 

5Q n = [ ( A T E A+yl)- 1 A T E' 1 (5R+ A8Q) ] n _ x 

SRn = R m -R c (Q n ) -> iterate to minimum. 

5Qn = Qn-Qo (but 8 t = O.lT, 5P cld = P cld (k)-P cld (k-1)] 


n = iteration number; k = pressure (P) grid number. 

Q = variables [i.e., T(P), q(P), Ts, Ps, Pcld, t, De(x) ...] 
Q 0 = regression retrieval. 

= measured radiance 
R C (Q) = calculated radiance from Q 
A = Jacobian Matrix 

Ay = [R i (Qj)-R i (Qj-8Qj)]/ 5Qj (for cloud parameters) 
E = error covariance matrix 

y = Lagrangian multiplier (optimized in the iteration) 


Simultaneous 1-D Var. Iterative Matrix Inversion 
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Radiative Transfer Model (include clouds) 
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= upwelling spectral radiance at the top of atmosphere 
= cloud transmissive function 
= cloud reflective function 
= upwelling emission below the cloud 
= downwelling emission below the cloud 
= emission from the cloud 
= upwelling emission above the cloud 
= downwelling emission above the cloud 

= surface emissivity 
= Planck function 

= total transmittance from any given level to an upper 
boundary such as cloud level or the top of the atmosphere 
= the total transmittance from any given level to a lower 
boundary such as cloud level or the Earth’s surface 
= transmittance between the cloud level and the Earth’s 
surface 

= transmittance between the top of the atmosphere and 
cloud level 
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EOF Physical Regression Inversion 


Statistics are formulated for one class of data 
which contains all cloud height conditions 

and 

2 other classes for which the cloud phase has 
been stratified to liquid and ice. 
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R = radiance 

P s = surface pressure 

S = number of sample profiles 

= radiance deviation from the mean 
M = covariance matrix of 
E = eigenvectors of M - EOFs 
C = radiance EOF amplitudes 
A = {T s , \|/, T, q, ... Held, xcld, De, Pha} parameters 
K = regression coefficients 
\| / = emissivity EOF amplitudes 
8 = emissivity 
e = emissivity eigenvectors 

F x = cloud transmissive function {Held, xcld, De, Pha} 

F r = cloud reflective function {Held, xcld, De, Pha} 

Held = cloud hight 

xcld = cloud optical depth 

De = cloud particle diameter 

Pha = cloud phase (ice or water cloud) 
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Global Training for LaRC Algorithm 
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Channel Used in LaRC Retrieval Algorithm 


IASI: 5008 channels for regression, 1697 channels for physical retrieval 
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Retrieval Parameters from this System 


Brightness Temperature or Radiance Spectrum 


T: temp. prof. 
Surf: surf. prop. 
Cld: cloud prop. 
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Retrievals under clear conditions: 

• (^ urfac ^)properties (skin temp and 
emissivity). 

• Atmospheric temperature and<j noistur e) 
profiles. 

• Atmospheric CO and 0 3 abundances. 


Retrievals under cloudy conditions: 

• Atmospheric profile through optically thin 
cirrus clouds and above optically thick 
clouds. 

• Effective( clou3 )parameters (i.e., cloud top 
pressure, particle size, and optical depth). 
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***** IASI Retrieval Demo: Cloud Top Height 
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IASI vs. GOES-12: Cloud 
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IASI Retrieval Demo: Cloud Optical Depth 


Cloud Optical Depth 
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IASI Retrieval Demo: Cloud Particle Size 


Cloud Particle Diameter (^m) 
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IASI Retrieval Demo: Surface Skin Temp 


Surf Temp (K) 



-150 -100 -50 O 50 100 150 

Longitude (deg.) 


250 260 270 280 290 300 310 320 


Dan Zhou @ ITSC-16; May 7, 2008 




Langley Research Center 


N ^®f IASI Retrieval Demo: Moisture Distribution 


HurricaneN RH Horizontal Distribution (%) 24.49 km 
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Hurricane Dean (2007.8.20) 




This image of Hurricane Dean is a 
combination of observations from 
NASA and NOAA satellites. The 
clouds were observed by a NOAA 
GOES-12 at 20:45 UTC on 20 th 
August 2007. The land surface is a 
summertime image from the NASA 
Blue Marble image collection. 


Dan Zhou @ ITSC-16; May 7, 2008 



. iLli'i.drc Jikg.i 


Langley Research Center 


N'VSA 


Hurricane Dean Observed with IASI 
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Hurricane Dean Observed with AIRS 
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IASI vs. AIRS: Cloud Top Height 


Metop Satellite / IASI 
Time: -15:14 UTC 
Eye Location: 17.8 N, 82.5 W 
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Aqua Satellite / AIRS 
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Joint Airborne IASI Validation Exp. 


Location/dates: 

Ellington Field (EFD), Houston, TX, 14 Apr - 4 May, 2007. 
Aircraft: 

NASA WB-57 (NAST-I, NAST-M, S-HIS); 

UK FA AM BAel46-301 (ARIES, MARSS, SWS; dropsondes; in-situ 
cloud phys. & trace species; etc.). 

Satellites: 

Metop (IASI, AMSU, MHS, AVHRR, HIRS). 

A-train (Aqua AIRS, AMSU, HSB, MODIS; Aura TES; CloudSat; 
and Calipso). 

Ground-sites: 

DOE ARM CART ground site (radiosondes, lidar, etc.) 
Participants: 

include NASA, UW, MIT, IPO, NO A A, UKMO, EUMETSAT, 
ECMWF, ... 


UK FAAM BAe 146-300 


NASA WB-57 
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Case Study and Validation (2007.04.29) 
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IASI Retrieval: Cloud & Surface 



Cloud Parameters: 
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IASI Cloud vs. GOES Image (4.29.2007) 
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IASI Retrieval: ATemp and RH Fields 
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IASI Retrievals vs. Radiosondes 
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IASI Regression vs. Physical Retrieval 


Statistical Regression Retrieval: RH 
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1. The retrieval improvement based on the EOF statistical regression through 
physical iterative retrieval is only contributed by IASI measurements as the 
minimum information methodology used. 


2. A high-vertically-resolved atmospheric structure is captured very well by IASI 
measurements and/or retrievals; not only in the troposphere, but also in the 
boundary layer. 
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IASI (15:48 UTC) vs. AIRS (19:30 UTC) 
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Summary and Future Work 


1. A state-of-the-art IR-only retrieval algorithm has been developed with an all- 
seasonal-global EOF Physical Regression and followed by 1-D Var. Physical 
Iterative Retrieval for IASI, AIRS, and NAST-I. 

2. The benefits of this retrieval are to produce atmospheric structure with a 
single FOV horizontal resolution (~15 km for IASI and AIRS), accurate 
profiles above the cloud (at least) or down to the surface, surface parameters, 
and/or cloud microphysical parameters. 

3. Initial case study and validation indicates that surface, cloud, and atmospheric 
structure (include TBL) are well captured by IASI and AIRS measurements. 
Coincident dropsondes during the IASI and AIRS overpasses are used to 
validate atmospheric conditions, and accurate retrievals are obtained with an 
expected vertical resolution. 

4. JAIVEx has provided the data needed to validated retrieval algorithm and its 
products which allows us to assess the instrument ability and/or performance. 

5. Retrievals with global coverage are under investigation for detailed retrieval 
assessment. It is greatly desired that these products be used for testing the 
impact on Atmospheric Data Assimilation and/or Numerical Weather 
Prediction. 
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IASI vs. AIRS: Skin Temp. 
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Eye Location: 17.8 N, 82.5 W 
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